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Abstract 

We extend the investigation of BPS saturated t’Hooft-Polyakov monopoles in x to the 
general case of SU{N) gauge symmetry. This geometry causes the resulting N — 1 coupled non¬ 
linear ordinary differential equations for the iV — 1 monopole profiles to become autonomous. One 
can also define a flat limit in which the curvature of the background metric is arbitrarily small 
but the simplifications brought in by the geometry remain. We prove analytically that non-trivial 
solutions in which the profiles are not proportional can be found. Moreover, we construct numerical 
solutions for N = 2,3 and 4. The presence of the parameter N allows one to take a smooth large N 
limit which greatly simplifies the treatment of the infinite number of profile function equations. We 
show that, in this limit, the system of infinitely many coupled ordinary differential equations for 
the monopole profiles reduces to a single two-dimensional non-linear partial differential equation. 

1 Introduction 

The subject of’t Hooft-Polyakov monopoles DP m has been extensively studied in the literature, (see 
for reviews m m 15]). The fundamental role these objects are believed to play in modern theoreti¬ 
cal physics ranges between all branches of the subject (for a comprehensive review see, for example, 
my These solitonic objects have become predominant characters of most modern-day non-Abelian 
gauge theories, including supersymmetry and string theory (see for example [8|) despite having 
never been seen in experiments. One aspect of their description which is less studied is their intimate 
link with the topology of the underlying space. This paper is an extension of our previous work [9| 
which investigates solitonic t’Hooft-Polyakov monopoles in a cylindrical topology with spherical cross- 
sections (for related work see also |10j |lljl. In our previous investigation we showed that for SU{2) 
monopoles living in this topology all field profile equations become autonomous (that is, they don’t 
depend on the radial variable r explicitly) which allows one to treat them by analytical tools normally 
not available. In this work we extend this set-up to the general case of SU{N) gauge symmetry. The 
main mathematical tool used to achieve this is the formalism of harmonic maps which is reviewed 
in section 2. We will show that, as usually happens, by introducing an arbitrary dependence on the 
parameter N, one can find find a simplifying limit in which N becomes large. The simplifications 
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brought about by this limit are especially treatable in the harmonic map formalism and constitute the 
main reason for choosing this over other methods. 

The paper is organized as follows: section 2 is devoted to introducing the system with its La- 
grangian and the aforementioned topology. This section also includes a review of the harmonic map 
formalism used to describe SU{N) monopoles. In section 3 we discuss in detail the consequences of 
the large N limit and in section 4 we find solutions numerically for the cases N = 2,3,4. Finally we 
provide some conclusions in section 5. In the appendix we provide a rigorous mathematical proof that 
the resulting equations for field profiles always allow non-trivial solutions (given a suitable bound on 
the shape of the cylinder). 


2 The System 

We begin this section by reviewing the general algorithm to obtain spherically symmetric SU{N) 
monopole solutions using rational maps of the Riemann sphere into flag manifolds. We refer the 
reader to [12], on which this review is based and from which we borrow our notation, for the relevant 
mathematical details. 


The action of the SU(N) Yang-Mills-Higgs system in four dimensional space-time is 

5ymh = I Tr ^ , (i) 

where the Planck constant, the speed of light and the gauge coupling constant have been set to 1. 
The remaining dimensionless coupling constant is A. In this equation, is an su{N) valued scalar 
field, the covariant derivative is Di = di + [Ai,] and Fij = diAj — djAi 4- The SU{N) BPS 

monopoles are finite energy solutions of the BPS equation 

A'h = (2) 

obtained by an energy minimization argument from eq.([^ when A = 0. In polar coordinates with 
z = e®‘^tan(0/2) the flat space-time metric reduces to 

ds'^ = —dt'^ -|- dr^ + 7 - 1 .„.„ dzdz, (3) 

(1 -|- \z\^Y 

and one can use the ansatz 


^ = -iAr = A^ = H-^d^H, Ag = 0, 

for H G SL{N,C), to reduce the matrix system of BPS equations to 

| 2\2 


(l + \z\^) 

dr {H-^drH) + ^= 0 . 


(4) 


(5) 
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The ansatz used in eq.Q (which is the most natural generalization of the ’t Hooft-Polyakov 
hedgehog ansatz for SU{2)) describes three-dimensional topological defects with non-Abelian magnetic 
charge given by 

Qm= [ (6) 

Jsr 

where Sr denotes a sphere centered around the monopole and Bf = \e-ijkFjk is the non-abelian mag¬ 
netic field. More generally, the dimensionality of a topological defect (see, for a detailed explanation, 
m and [I3]) is determined by the homotopy class of the corresponding ansatz. In general, for systems 
living in a spacetime A4 and with a gauge symmetry G broken down to a subgroup H one can label 
the toplogical charges of solitonic solutions by the degree of maps of 5” G Ad into the coset space 
GIH, 

7^n{G/H), (7) 

where n varies according to A4 and which pattern of symmetry breaking is considerecQ For the case 
of SU{N) monopoles considered here, we will be interested in the homotopy class 

TTS {SU{N)/J) , (8) 

where J can vary between U{N — 1) and U (1)-^“^ corresponding to how the gauge symmetry is broken. 
The case U{N — 1) corresponds to minimal symmetry breaking whilst 17(1)^“^ is the case of max¬ 
imum symmetry breaking and these depend on the vacuum expectation value of <h. These solitonic 
objects are stable unrelated solutions of different energy minimization requirements, however one of 
the interesting outcomes of the present analysis is that different topological objects can actually be 
difficult to tell apart (at least by looking at the equations of motion) when such topological objects 
are analyzed within space-like regions with non-trivial topology. In particular, we will show that the 
field equations for non-Abelian BPS monopoles (possessing non-trivial second homotopy class as per 
eq.Q) within the bounded tube-shaped region defined in the next section (see M) are related to the 
field equations of domain-wall objects in flat topology and with a similar homotopy class eq.Q (see 
|14j l through a simple field redefinition. 


In order to find spherically symmetric monopole solutions to equation Q one can follow a simple 
algorithm: first, one needs the spherically symmetric maps into (see |12] ) which are given by 


fj=Z^ 





(9) 


where the expression in the square root denotes the standard binomial coefficient. Then using the A 

^For instance, in flat space, vortex-like objects have a non-trivial first homotopy class 7 ri([/(l)) , monopole-like objects 
(which are point-like in space) have a non-trivial second homotopy class 'K2{SU{2)/U{!)) while instanton-like solutions 
(which are point-like in space-time) have a non-trivial third homotopy class tt3{SU{N)). 
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operator defined as 


A/ = dj 


( 10 ) 


t |^|2 , 

and applying it iteratively A^f = A (A^“^/), A: = 0,...,A — 1, one can construct a projector matrix 


If 

I/P 


Pk = P{Af), P{f) = ff, (11) 


satisfying P^ = P = P^ which is used to parametrize the general N x N SL{N, C) Hermitian matrix 
H appearing in eq.([^ by 

where gi are general profile functions which depend only on r and the 1/N factor denotes the identity 
matrix divided by N. Substituting eq.(12) into eq.([^ gives a general matrix of equations which can 
be decoupled for each gi (here' denotes differentiation w.r.t r). More generally, following [T^ one can 
find a convenient form of the resulting equations in terms of N and I = 0,N — 2, the index labelling 
the profile function gi, 

^ ^ ^(i + l)Fi + 2j]^Fi-^(/ + l)(A-^-l)(exp(FO-l) = 0, (13) 


N 


i=0 


fc=0 i=k 


where Fi = gi — gi+i and Fn -2 = 9N- 2 - Decoupling these equations gives N — 1 equations for the 
profile functions, we include below the examples for SU{N) with N = 2,3, which are respectively 


(7o + ^(l-e9«) = 0, 


(14) 


and 


- 91-I («”-!) += 


(15) 


(16) 


As shown in m, these equations have analytic solutions corresponding to magnetic monopoles. 
In the simplest case of A = 2, the solution is 


go = 2 log (2r/ sinh 2r), 


(17) 


which describes a single BPS saturated monopole, with energy equal to its charge Qm- In the case 
where N > 2, the equations still allow for at least one analytic solution in which the profile functions 
are chosen proportional to each other. Indeed if, for example, we choose N = 3 then there exists the 
solution with go = 2gi. In general, one has a solution for all field profiles proportional when 


9i 


N-{i + l) 


= 9 


(18) 
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with g given by eq. (© and i = 0 ,N — 2. 

2.1 SU{N) monopoles on x 8“^ 

In this section we discuss a simple geometrical modification of the above set up and its consequences. 
We wish to consider the above system in x (or Al^ x 5^ x with metric 

ds^ = —dt^ + + R^{d6‘^ + sin^ 9 dcf^) , 0 < r < L , (19) 


O<0<7r, 0<(/><27r, 


( 20 ) 


where L is a longitudinal length and i?o is a constant with the dimension of length related to the 
size of the transverse sections of this topology. In complexified coordinates, those appropriate to this 
paper, the metric reads 


ds^ = —dr + dr + 


ml 


^dzdz. 


( 21 ) 


(l + |z|2)2 

This metric describes a tubular geometry with spherical caps as cross sectional slices. The non- 
vanishing components of the Riemann tensor of this space are proportional to I/Rq'- 


Rpi/pa ~ 


1 


( 22 ) 


This simple modification of the geometry leads to a dramatic simplification in the resulting equa¬ 
tions for the profile functions. As shown in [9], for the case of SU{2), the resulting energy minimization 
equations for the monopole profile functions are autonomous, that is they don’t involve any explicit 
powers of r. Below we find that this also happens for the SU{N) case. As shown later in the paper, 
one can easily define a flat limit in which the curvature of the background metric is as small as one 
wants keeping, at the same time, all the simplifications brought in by the background metric in eq. 


(21) (two different ways to achieve the flat limit will be described in the following sections). Moreover, 
the present formalism introduces the possibility to define a large N limit. We will consider this in 
detail in a further section. 

Equation when analyzed in the background metric given by eq.(|21[) is modified to 


dr {H-^drH) + 


Rl 


d, {H-^d,H) = 0, 


(23) 


and equation (13) becomes. 


N-2 


I N-2 




2 (Z + 1) 


i=0 


k=0 i=k 


Rl 


(Ar_;_i)(exp(Fz)-l) = 0, 


(24) 


where, as in the previous section, Fi = gi — The equations now become autonomous. In line 


with the discussion around eq.(18) we find that for N > 2 there exist analytical solutions in which all 


the profile functions are proportional. However, in the following we will show both numerically and 
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analytically that non-trivial solutions in which the profiles are not equal (and hence do not correspond 
to trivial embeddings of SU{2) into SU{N)) also exist. 


Let us take a moment here to connect the previous results (and the more standard monopole 
notation) to the cnrrent notation. For the simplest case of SU{2) the equation reduces to 


<?o + ;^(l-e"°) = 0. (25) 

It was shown in |9] that for N = 2 nsing standard polar coordinates (r, 9, (j)) and an ansatz of the form 

Af, = ikir)-l)U-^d^U , ^ = i>ir)U, (26) 

U = nHi , U-^ = U , (27) 

= sin 9 cos cj) , = sin 9 sin cj) , = cos 9 , (28) 

where the f are the standard Panli matrices, the BPS equations reduced to 

dr^l^ + = 0, (29) 

Kq 

drk — k-ip = 0, (30) 

for which a general solution was proposed of the form 


k = exp(u) , 


(31) 


= drU , 

where the function u{r) is the inverse of the following integral 


pu{r) 

J 1 /( 0 ) 


2 ( Iq + 


exp {2z) — z 


- 1/2 


dz = ±r , 


(32) 


(33) 


with Iq an integration constant. This is consistent with the above construction, as one indeed expects. 


upon the identification u = qq!2 where eqs.(29) and (30) reduce to eq.(25). 


The autonomous set of equations (24), of which eq.(25) is both the N = 2 and the degenerate equa¬ 
tion for proportional fields representative, has some interesting properties which we discuss throughout 
the paper. One of which is that it coincides with the equation for a domain wall separating the Higgs 
and Coulomb phases in the Abelian - Higgs modeQ In [T4] , the equation for the scalar field describing 
the domain wall was found to be 

(log(/>)"^ = (())2-n^) , (34) 


^We thank S. Bolognesi for pointing this out to us. 
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where e is the U{1) gauge coupling and is the constant appearing in the quartic potential (similar to 
that appearing in eq.([^). The above equation is a particular case of the well known Taubes equation 


m- Upon identifying cj) = v exp{gQ/2) and = 1/Rq we see that eq.(|34[) becomes eq.(|25[). What 
is intriguing about this observation is that topological objects which are quite different (possessing 
different non-trivial homotopy classes) can be described by identical solutions if they are constrained 
to live in space-time regions with non-trivial topology. In fact, since there is always a solution of the 


N equations in which the field profiles are proportional as in eq.(18), this relation can be trivially 
extended to monopoles of topological charge N. 


3 The large N limit 

The large N limit introduced in Yang-Mills theory in m (see also |20] and [2T]; for two detailed 
reviews see [22]) is a very powerful tool to analyze non-perturbative features in gauge theories (such 
as confinement, bound states and so on). The non-trivial scaling behavior of physical quantities within 
the large N expansion arises from the fact that Feynman diagrams have different weights depending 
on the topology of the surfaces they can be drawn on (once N is considered as a large number). Thus, 
this non-trivial behavior is purely “quantum” in nature and one would expect that no non-trivial large 
N behavior should be found when analyzing classical BPS equations (as we are doing in the present 
paper). In fact, the present formalism shows that a non-trivial scaling with N emerges already at the 


level of the BPS field equations eq.(24). Naively taking N very large means increasing the number 


of non-linear coupled equations, or equivalently the number of monopole profiles, which seems like an 
un-neccesary complication. 

Fortunately, there is another dramatic simplification in this limit which makes this problem treat¬ 
able. The key issue is that in this limit, through an appropriate ansatz of the interpolating discrete 


function, one can replace the discrete label I = 0,..., — 2 in eq.(24) by a continuous label. The main 

reason is that in the limit in which N —)• oo, where the range of values of I becomes infinite, each 
discrete jump in its value (which is just 1 of course) becomes infinitesimal compared to the range. In 
other words, were we to rescale each step A/ = 1 of the complete range by N, this would become 
infinitesimal in this limit. 


Let us consider this simplification in detail. First we take eq.(24) and replace the discrete labels 


Fi by a discrete function F{i) (we omit the dependence on the radial variable for simplicity). This 
is simply a renaming which is useful to understand the subsequent manipulations. Clearly, since i is 


integer, this is a function on Z. Then we multiply eq.(24) by 1/N^ and define the quantity 


X = l/N. 


(35) 


Then, in the limit in which Y —>• oo, the discrete label I can be replaced by the continuous variable 
X G Consequently we can replace the discrete function F{1) living in Z by the continuous 
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function F(x) living in The equation then becomes 


v) g V E v) g) -'» 


N ^ N 

k=0 i=k 


m 


Now let us manipulate the sums, define 


(36) 


i/N = y, z = k/N, p = i/N, 


(37) 


then simple manipulations lead to 

l-l/N 

g 


X ^ 1-2/N ^ 


d*+v) E (!/+l/W)iF(JV9) + 2X;^ ^ XPlNp) 

y=0 z=0 p=z 


- J 2 + ^) (1 - (^ + ^)) (exp(F(lVx) - 1)) = 0. 

Finally, taking the large N ^ 00 limit, in which we can safely replace the sums by integrals 

l-l/N 

I 7 [G(!/)I 

lO 


L-V^ . .1 

[Giy)] / dy G{y) , 

^^ iV N^oo /n 

y =0 -^0 


(38) 


since, from eq. (35), 1/N plays the role of “dx” in the mathematical definition of Riemann and Lebesgue 
integrals, expression eq.pSl) greatly simplifies and becomes 


-2x y 
Jo 


52 

^G{y,r) 


X ri r q2 


dy + 2 r [ 

Jo Jz L 


^G{y,r) 


2x 

dydz - -^{1-x) (exp(G(x, r)) - 1) = 0, (39) 

Kq 


where we defined G(x, r) = F{Nx, r): G{x, r) will be denoted as the complete profile function since it 
encodes information about all the elementary profiles Fi{r) at the same time. 

From here throughout the rest of the paper (unless specified) we switch to dimensionless units 


p = vr, Rq = vRq, E = E/v, 


(40) 


where v is the parameter appearing in the potential with dimensions of mass and E is the energy of 


the solution. Then eq.(39) becomes 


-2x y 
Jo 


92 


X /•! ■ 


dy + 2 
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2x 

dydz - 2(1 “ ^) (exp(G(x, p)) — 1) = 0 . (41) 

Ro 


Taking two derivatives with respect to x gives single non-linear elliptic partial differential equation for 













the complete profile function G{x, p) which reads 


dp^ dx‘^ 


x{l — x) 
Ro" 


(exp(G(x,p)) - 1) 


= 0 


(42) 


Solving this equation with the boundary conditions equal to the profile functions in the radial direction 
is expected to yield a function which interpolates smoothly between the monopole prohles. This is a 
non-trivial numerical task which we defer to a later publication. The main important advantage of 


eq. (42) over the original system in eq.(24) is that instead of having, in the large N limit, infinitely 


many coupled equations for infinitely many unknowns one has just one non-linear partial differential 


equation for the complete profile function G{x,p). Eq. (42) is well suited for numerical analysis and 
one could also easily apply the techniques described in the previous section to prove the existence of 
solutions with suitable properties 

It is worth emphasizing here that the only assumption used to derive eq.(|39|) starting from the 


original system in eq.(24) is that there exists a smooth interpolating function F{x) in which the group 


theoretical label ^ becomes a continuous variable x. A necessary condition is the existence of a large 
N limit of the field theory. Although the proof of this result is not available yet, there is a huge 
amount of evidence in the literature (see j22| and references therein) that it is actually possible to 
define a smooth large N limit in many different areas (from gauge theories to matrix model and so 
on). Whether the profile functions Ffir) as functions of the discrete label I follow a pattern regular 
enough to be interpolated by a continuous (smooth) function depends on the gi solutions at large N. 
Therefore, as used in this paper, this assumption translates to an ansatz on the solutions we consider 
in the large N limit which, in detail, involves the function Ffir) to be doubly differentiable. This 


ansatz leads directly to eq.(42). Remarkably the above argument can also be directly extended to the 


case in which one starts off with a flat topology and considers standard SU{N) BPS monopoles in flat 
space. 

We will now discuss two different ways to consider the flat limit in the background geometry in 
I ~ 2 

eq.(21) in which Rq —)■ oo. 


3.1 The flat limit as a perturbative limit 


Another important advantage of eq.(42) over the original system eq.(24) is that, in the large N limit, 

" —2 

it discloses the role of the adimensional curvature parameter Rq as coupling constant of the master 


equation given by eq.(42) in such a way that the flat limit corresponds to a “weak field” limit. 


Indeed, from eq.(42) it is clear that in the large Rq limit (in which the curvature vanishes) the 

" 2 

equation for G (x, p) becomes just a linear partial differential equation while the limit in which Rq is 

small (so that the curvature is large) the non-linear effects become strong as well. Thus, this formalism 

provides one with a clear perturbative scheme in which the equation for the complete profile function 

~ 2 

G {x,p) becomes linear in the flat limit in such a way that the curvature parameter I/Rq plays the 
role of a coupling constant. 
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In fact, if one would put the monopoles on a flat geometry from the very beginning the corre¬ 


sponding system would be eq.(13) and such a system does not admit any obvious perturbative scheme 


in the large N limit although one would arrive at an equation similar to eq.(42). Indeed, following 


similar steps to the previous section, if one would start from eq.(13) assuming in the large N limit a 


continuous dependence on the group label I, then one would arrive at the following large N limit for 
the complete profile function Gf{x,p): 


dp^ ^ dx'^ 


x{l — x) 


(exp(GF(x,/9)) - 1) 


= 0 , 


(43) 


where the label F has been added to emphasize that Gf{x,p) is the complete profile function of the 


system of monopoles described by eq.(13) which live, from the very beginning, on a flat metric. The 
difference between eq.M2| and eq.(43) is then apparent: in the former equation Rq ^ plays the role 


of coupling constant allowing a perturbative analysis of the equation while in the latter equation one 
cannot do this since Rq has been replaced by p~‘^ where p is one of the independent variables of the 
equation. 


3.2 The flat limit as a geometrical bound 

An alternative way to define a flat limit for the background metric corresponds to rescaling the original 


longitudinal variable in eq.(24) such that 


p = r/Ro . 


(44) 


Since the length of the tube-shaped region in which the monopoles are living is L, the above rescaling 
is equivalent to considering a tube of adimensional length 




(45) 


Then, one can take the flat limit taking Rq and L simultaneously large in such a way that L stays 
finite: 

Rq ^ oo , L ^ oo \ L = const / 0, oo . (46) 


In this limit, the master equation eq.(42) for the complete profile function simply becomes 


^ ^ (exp(G(x, p)) - 1)] = 0 . 


(47) 


As shown in the appendix it is important to note however that, in order to be sure that when con¬ 
sidering the flat limit in eqs.(|44[), ( |45[ ) and ( |46[ ) non-trivial solutions always exist, one should take the 
limit in such a way that the inequalitie^ in eqs.(85) and (86) are never violated. This fact can be 
interpreted as a sort of bound on the shape of the cylinder which cannot be too “slim” since L ^ 1 


^In the generic SU{N) case, the application of the Schauder theorem would give very similar inequalities. 
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Figure 1: BPS monopole profile function for SU{2). The plot corresponds to the numerical value (in 
appropriate units of u) i?o = 1-45. 


would violate eqs.(85) and (86). 


4 Numerical Solutions 

In this section we provide numerical solutions for the profile functions in the specihc cases of = 2,3 
and 4. For every value of Rq there exist a unique solution for the profile functions with an integer 
value of Qm- The numerical strategy is therefore the following: we provide boundary conditions for 
the profile functions and vary Rq until a solution with the desired topological charge is found. This 
guarantees that the solution obtained is a solution (for a given Rq) which represents a BPS monopole 
with topological charge Qm- The finite length cylindrical topology is implemented numerically by 
imposing the boundary conditions at a finite cutoff. The numerical procedure is a second-order finite 
difference procedure with accuracy 0(10“^). 


4.1 N = 2 


Let us proceed to solve equation (25) numerically. We fix uL = 50. We fix the boundary conditions 
on the profile function to be 

5o(0) = 1, ffo(50) = 0, (48) 


and we find that at i?o = 1-45 the energy reads 


E = 


Qm 

Att 


(e9o _ 1) ^ 
^ ’ 2 


p=50 

p=0 


1 . 000 , 


(49) 


which is the expected one-monopole solution. This solution is shown in figure 1. 

For this value of Rq one can use the topology of the system to find the solution corresponding 
to placing the monopole at the other end of the cylinder. This solution is found by imposing the 


11 












goip) 

1.0 

0.8 

0.6 

0.4 

0.2 



30 40 



P 


Figure 2: BPS monopole profile function for SU{2) with boundary conditions at other end of the 
cylinder. The plot corresponds to the numerical value (in appropriate units of v) Rq = 1.45. 


conditions 

ffo(O) = 0, 5 o(50) = 1, (50) 

and is shown in figure 2. For this solution, which is just an inversion of the previous solution about 
the center of the cylinder, the energy is the same, as expected. Note that this is not an anti-monopole 
solution as the topological charge is the same. 


At this point a comment is in order regarding the SU{2) solutions of [9]. Below eq.(33) we pointed 
out that there is a direct map between the profile function go and the solutions obtained without 
using the harmonic map formalism. Namely we showed that using u = go/2 one can map the two 
first order BPS equations into the profile equation for go- However, the solution for go presented in 
figure 1, when translated to the standard Higgs and gauge fields i/ = 5o/2 and k = exp{go/2) does 
not reproduce the solution found in [9] even though both solutions have the same topological charge 
(the reader may argue that the values of Ro used for both plots are not the same, but this comment 
applies also when these values are made equal). This is however expected, in order for both solutions 
to match, since A:(0) = 0, one would have to provide a singular boundary condition for the profile 
function g'o(O) = —oo, which is numerically impossible. The non-uniqueness of the solution with 
given charge is made possible because in this topology one does not require that the prohle function 
be regular at the radial origin, i.e. = 0 precisely because (unlike what happens in flat metric in 
spherical coordinates) there is no preferred origin. 
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4.2 N = 3 


In this section we wish to solve 


- 5 ^' + ^ 


9i 


2 

2 


Ro 


1) + ^ - 1) = 0 

Ro 

(51) 

1) + ^2 - 1) = 0. 

Ro 

(52) 


where ' denote differentiation w.r.t p. In Figure 2 we show the solution corresponding to the case where 
the profile functions are proportional. This is easily obtained by demanding boundary conditions of 
the form 

5o(0) = l, ffi(0) = 0.5 (53) 


50 ( 50 ) = 51 ( 50 ) = 0 . 


The general energy equation reads 

Qm 


E = 


dvr 




p=50 


p=0 


(54) 


(55) 


For the solution in Figure 3, by construction we find numerically 

E = 1.0006. 


(56) 


We can however also look for solutions which are not proportional. These are more interesting solutions 
from the numerical point of view as they are full solutions of the coupled equations rather than the 
reduction of all of these to one equation for a single profile function. 

Using the topology of the system we may also find solutions which have non-vanishing boundary 
conditions at both ends of the cylinder and therefore not proportional. In this case we seek solutions 
with boundary conditions of the form 


5o(0) = l, 5'i(0) = 0, 
5(,(50) = 0, 5 i(50) = 1, 


(57) 

(58) 


these are shown in Figure 4 for a very specific choice of Rq. The energy for this solution, calculated 

E = 2.001, (59) 


numerically using eq.(55) is 


and hence it is tempting to identify this solution as a monopole-anti-monopole state, with each species 
of monopole living at each end of the cylinder. However once again the boundary conditions lead to 
monopole charges which are not consistent with this picture. 
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Figure 3: Proportional monopole profile functions for SU{3). The plots correspond to the numerical 
value i?o = 1- For these plots go = 2gi. go is the solid line whilst gi is the dashed line. 




P 


Figure 4: A solution with non-vanishing boundary conditions at both ends of the cylinder at Ro = 4.28. 
Dashed line is gi, solid line is go- 
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Figure 5: An S?7(4) two monopole solution at Ro = 5.9. Thin dashed line is g 2 , medium dashed line 
is gi and solid line is i^o- 


4.3 N = A 

In this case we wish to solve the system of equations 


-9o + ^ - 1 ) + ^ - 1 ) = 0 , 

Ro Ro 


9i 


Ro 


^ (e9o-9i - 1) + A - 1) + ^ - 1) = 0, 


Rq 


Ro 


- 92 - 4 ^ - 1 ) + ^ - 1 ) = 0 , 

Ro Ro 


(60) 

(61) 

(62) 


and the energy evaluates to 
Qm 


E = 


dvr 


= [3e5o+^2 ) + 4^2,1 ^ (_3^/ - + (1 + 3e^^) )] 


p=50 


p=0 


(63) 

We look for the analogous two-monopole solution of the previous section. Therefore we impose 
the boundary conditions 

5o(0) = l, ff'i(0) = 0, 52(0) = 0, (64) 


^(,(50) = 0, 5i(50) = 1, 5^(50) = 0. 


(65) 


The corresponding solution shown in Figure 5 is found to have E = 2.000 at Rq = 5.9. If one calcu¬ 
lates the topological charge at each end of the cylinder one once again does not find equal contributions. 


In each case, it appears that the profile functions concentrate around the end-points of the cylin¬ 
drical topology and want to vanish in the intermediate region. 
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5 Conclusions 


In the present paper, using the harmonic map formalism, we have extended the investigation of BPS 
saturated t’Hooft-Polyakov monopoles in x S'^ to the general case of SU{N) gauge symmetry. We 
found that, as per the N = 2 case investigated previously in [9], all equations for the monopole profile 
functions become autonomous. For some specific cases we solved these equations numerically and in 
all cases we demonstrated analytically the existence of non-trivial solutions in which the field profiles 
are not proportional. Furthermore, we investigated the remarkable power of the large N limit in this 
scenario, where we have shown that one can, under a suitable condition of smoothness, reduce the 
infinite set of profile equations to a single partial differential equation for an interpolating function. 
This equation encodes the subtle role played by the curvature parameter Rq and elucidates the flat 
i?o oo limit both from a perturbative and a geometrical aspect. We leave the numerical treatment 
of this equation to further work. 
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A Existence of solutions 


In this appendix we will provide a rigorous mathematical proof that solutions to the set of equations 


obtained from eq.(24) exist even in the case that the fields are not proportional. This procedure does 


not yield the solutions themselves, these are shown numerically in a previous section of the main text, 
but is nonetheless an important part of the analysis of the field profile equations and the space of their 
solutions. In particular, one of the inequalities derived in this section will be useful in the discussion 
of the flat limit. 


The basic mathematical tool used here is the Schauder theorem (see for a detailed pedagogical 
review m ). For simplicity, we will focus on the SU{3) case but the same argument can be easily 
extended to the general case. 

The statement of the Schauder theorem (HZl [T8] l is the following: let 5 be a complete metric 
(Banach) space so that a distance d {X, Y) between any pair of elements of the space is defined by 

d{X,Y)eR, X,YeS, (66) 

and such that, with respect to the chosen metric, from every Cauchy sequence one can extract a 
convergent subsequence (this is the inclusion of “complete” in the definition of the metric space). Let 
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C be a bounded closed convex set in S and let T be a compact operatoi]^ from the Banach space S 
into itself such that T maps C into itself: 


T[]-.C^C. (67) 

Then the map T [.] has (at least) one fixed point in C. In other words, under the above hypothesis, 
there always exist a solution to the equation 


T[X]=X . (68a) 

Our task is to determine under what precise conditions this applies to the system of equations 


obtained by decoupling the field profiles in eq.(24) for the case of = 3. To do so, let us rewrite 


the system in eqs.(15) and (16) (with the factor r replaced by i?o everywhere) as coupled integral 
equations: 


go (x) = ao + box + 

Jo Jo 

gi (x) = ai + bix + 

Jo Jo 


2 2 

{exp [go (p) - gi {p)] - 1} + [exp {gi {p)) - 1] 

Ro Ro 


4 2 

[exp {gi (p)) - 1] -(exp [go (p) - gi (p)] - 1} 

Ro Ro 


dpds , (69) 

dpds , (70) 


where a* and bi represent the initial data for the two profiles go and gi and their derivatives at x = 0. It 
is a trivial computation to show that the above system of integral equations is equivalent to the system 


in eq.(24) with N = 3. The system of eqs.(69) and (70) can be written as a fixed point condition for 
the following vectorial operator ^ acting on pairs of continuous functions (x) = {go (x) ,gi (x)): 

f [0, L] X [0, L] [0, L] x [0, L] , 
t {x) = {go {x) , <71 (x)) G [0, L] X C" [0, L] , 


(71) 


2 2 

(exp [go {p) - gi (p)] - 1} + [exp {gi {p)) - 1] 

Ro Ro 

4 2 

[exp {gi {p}) - 1] - {exp [go {p) - gi (p)] - 1} 

Ro Ro 


f [go, gi] = f[t (^)] = {To {x) , Ti (x)) , 

where 

To (x) = ao + box + f f 
Jo Jo 

Ti (x) = ai + bix + f f 
Jo Jo 

The fixed-point condition is then simply 

t{x) = f[t{x)] 

operator T from a Banach space into itself (see, for a detailed discussion, m DS]) is called compact if and only 
if, for any bounded sequence {xn}, the sequence {T{xn)} has a convergent subsequence. 


dpds , 

(72) 

dpds . 

(73) 


(74) 
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where the operator ^ has been defined in eqs. (71), (72) and (73). Whilst this proves that a fixed 
point operator condition can be found, in order to apply Schauder’s theorem we have to show that ^ 
is a compact operator from a bounded closed convex sub-set of a Banach space into itself. This is not 
a trivial task, let us begin by defining the following metric in the space [0, L] x [0, L] (which is 
the Cartesian product of the space of the continuous function on [0,L] with itself): 


sup \go{x) - ho{x)\ + sup \gi{x) - hi{x)\ , 
^ x&[0,L] a;S[0,L] 

't = {9o{x) ,gi{x)) , 't = {ho{x) ,hi{x)) . 


(75) 


With respect to this norm, the space [0, L] x [0, L] is a Banach space (which we call S). 

Then we define a bounded closed convex sub-set C of the Banach space defined above (using the 
metric in eq.(75)) such that ^ maps C into itself, 

c ^{t {x) = {go (x) ,gi (x)) G S| Vx G [0,L] \go{x) - ao\ < B , \giix) - ai\ < B} , B , 

(76) 

where ai are the initial data appearing in eqs. (69) and (70) so that C is closed by definition. It is easy 
to see that C is bounded since 


•B > |5o(a;) - aol > |5o(a:)| - |ao| ^ |c/o(a:)| < B-h |ao| , (77) 

B > |5i(x) - ail > |5ti(x)| - |ai| =► |5'i(a^)l < B|ai| . (78) 

In order to prove that C is convex we have to check that if (x) and h (x) both belong to C then 

(x) -|- (1 — 0) h (x) also belongs to C G [0,1]. This is easily verified as 

jdgi (x) + (1-0) hi (x) - ail < 10 (gi (x) - ai)| |(1 - 0) (hi (x) - ai)| < 0B -k (1 - 0) B < 5 , 

(79) 

1050 (x) + (1-0) ho (x) - aol < |0 (50 (x) - ao)| -k |(1 - 0) (ho (x) - ao)| < 0B -k (1 - 0) B < 5 . 

(80) 


Now we can proceed to show that T is compact. 

First of all, we must show that if '^n (^c) = (gj} (x) (x)) is a sequence in C then the sequence 

(x)] = (To (n; x), Ti (n; x)) is uniformly bounded in C (namely, the absolute values of both 
components of ^ [ ~^ri (a^)] are bounded by a constant which does not depend on n hence ensuring 
that ^ (x)] belong to C Vn as well). Therefore we consider 


iTo (n;x)| = 


ao + box + 


rx PS 

Jo Jo 


Ti {exp W (p) - sf (p)] - 1) + T 2 [exp (pf (p)) - 1] 

flo Ro 


dpds 


< 
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dpds =► 


|®o| + 


\hL\ + r [ 

Jo Jo 


{exp bo" (p) - (p)] - 1} + ^ [exp (g^ (p)) - 1] 

Ro Ro 

2L2 


\To (n;x)| < jaoj + |6oL| + [{exp {2B + |ao| + jaij) - 1} + [exp (fi + |ai|) - 1]] , (81) 

Ro 


and, similarly, 


2L2 


[Ti {n;x)\ < jail + |6iL| + [2 [exp {B + |ai|) - 1] + exp [2B + |ao| + jaij] - 1] . (82) 

Ro 


In order to derive eqs.(81) and (82) we used that, because of eqs. (77) and (78) (which are equivalent 
to saying that the (x) belong to C V n), one has 


v?2 i9o (p) - 9i (p)] - 1} + v?2 [®^P i9i (p)) - 1] 
Ro Ro 


< 


Ro 


[{exp {2B + |ao| + jaij) - 1} + [exp {B + |ai|) - 1]] , 
Ro 

^ [exp (pr (p)) - 1] - ^ {exp bo" (p) - Pr (P)] - 1} < 
Ro Ro 

[2 [exp {B + bij) - 1] + exp [2B + |ao| + bij] - 1] , V n. 


(83) 


(84) 


Eqs. (81) and (82) show that, \i~in (x) = (po (x ), p" (x)) is a sequence in C, the sequence ^ (x)] 

is uniformly bounded. Moreover, one has to require that the sequence of images ^ [ ~^n (^c)] belongs 
to C as well. As always happens (see m and |18] ) this will give some constraints on the range on the 
parameters B, L and Rq. In order for the sequence of images to belong to C the following inequalities 
must be satisfied (as can be easily seen by comparing eqs. (@, ([^ and ([78|) with eqs. ([sT]) and 


(82)): 


2L? 

\boL\ + -^2 [{exp {2B + |ao| + |ai|) - 1} + [exp {B + |ai|) - 1]] < B , 
Ro 
2L^ 

\biL\ + [2 [exp {B + |ai|) - 1] + exp [2B + |ao| + |ai|] - 1] < B . 

Ro 


(85) 

( 86 ) 


These imply that in order for this theorem to work, the length L of the cylindrically-shaped region in 


which these non-Abelian BPS monopoles are living cannot exceed the bounds defined in eqs. (85) and 


(86). One cannot obtain a very large value for the allowed L by simply increasing B since the left hand 

I I I I ~ 2 

sides of eqs. (85) and (86) increase faster than the right hand sides, but the situation improves if Ro 


is very large, namely in the flat limit (in which case the exponentials containing B are suppressed). 
What is important however is that it is always possible the choose B, L and Rq in such a way that 


eqs. (85) and (86) are fulfilled. 


The next step to prove that T is compact is to show that if (x) = {go (x ), p" (x)) is a sequence 


19 

























in C then the sequence ^ [ ~^n (a^)] is equicontinuou^ To show this, we must evaluate, for a generic 
n, the absolute values of following differences: 


\To (n;x) - Tq (n;y)| = 


\Ti (n;x) - Ti (n;y)| = 


ry fs 

bo{x-y) + 

Jx Jo 

{exp K (p) - 9 f (p)] - 1) + -^ (exp ( 9 ? (p)) - 1] 

[Ro Ro J 

dpds 


(87) 

ry fs 

bi{x-y) + 

Jx Jo 

[e^P (p)) - 1] - v?2 {®^P bo ip) - 9i (p)] - 1} 

[72o Ro \ 

dpds 


( 88 ) 


where 0 < x < y < L. After some trivial manipulations (which use the fact that all the functions 
(x) belong to C and consequently eqs.(77) and (78) are satisfied) one arrives at 


\To (n; x) - Tq (n; y)\<\x-y\ 


\Ti {n;x) - Ti (n;y)| < \x - y\ 


2L 

l^ol + [{exp {2B + |ao| + |ai|) - 1} + [exp (fi + |ai|) - 1]] 

Ro 


2L 

\bi\ H-2 [2 [exp {B + |ai|) - 1] + exp [2B + |ao| + |ai|] - 1] 

Ro 


(89) 

(90) 


Thus, given any e > 0, we can choose 

5 < 


||bo| + |6i| + J4 [2 [exp {B + |oi|) - 1] + exp [2B + |ao| + \ai\] - 1]| 


(91) 


in such a way that both the choice of 6 in eq. . (91) does not depend on n and, 


\x - y\ < 5 ^ \Ti {n;x) - R (n; y)\ < - , V n, V i = 0,1 


(92) 


In summary, eqs. (81), (82), (85) and (86) show that, if ~i^n{x) = {gQ{x),gi{x)) is any se¬ 
quence in C, then the sequence ^ [ ~^n{x)] is uniformly bounded in C. Subsequently, eqs. (89), 
(90), (91) and (92) show that, if ~^ri (x) = (gQ (x) ,gi (x)) is any sequence in C, then the sequence 
T[tn (x)] is equicontinuous. Consequently, using the Ascoli-Arzela’ theorem (see [IZ])) from any 
sequence 7iy„ (x)] one can extract a convergent subsequence: this implies that the operator ^ is 
a compact operator from a bounded closed convex set into itself. 


Finally, the Schauder theorem ensures that eq.(74) (which is equivalent to our original system) 


has at least one solution. Moreover, it is always possible to choose appropriately the initial data a* 
and bi in such a way that the two profiles are not proportional. This concludes the proof on existence 
of solutions of the monopole profile equations. Moreover, one can also show by a similar procedure 
that the solutions are actually not just continuous but they also have continuous first and second 


®A sequence of functions {/n} is said to be equicontinuous if, given e > 0, 3 <5 > 0 such that \fn (x) — fn {y)\ < e 
whenever \x — y\ < 5 and, moreover, 5 does not depend on n (otherwise the sequence would be continuous but not 


equicontinuous: see EZ])- 
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derivatives|3 

The present rigorous argument can be easily extended to the SU{N) case with N > 3. Besides 
the intrinsic mathematical elegance of the fixed-point Schauder-type argument, the present procedure 
also discloses the presence of the bounds in eqs.(85) and (86) on the length L of the tube-shaped 
region in which these non-Abelian BPS monopoles are living. At the present stage of the analysis, it 
is not possible yet to say whether such a bound is just a limitation of the method or it signals some 
deeper physical limitation on the volume of the regions in which one constrains these non-Abelian 
BPS monopoles to live. Understanding whether or not such BPS monopoles can fit into very large 
cylindrically-shaped regions is certainly a very interesting question on which we hope to come back in 
a future investigation. 
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